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ABSTRACT 

We present results from neutral hydrogen (HI) observations of the tip of the Magellanic Stream (MS), 
obtained with the Arecibo telescope as a part of the on-going survey by the Consortium for Calactic 
studies with the Arecibo L-band Feed Array. We find four large-scale, coherent HI streams, extending 
continously over a length of 20 degrees, each stream possessing different morphology and velocity gra- 
die nts. The newly discovered streams provide strong support for the tidal model of the MS formation 
bv IConnors et alT ([20061 . which suggested a spatial and kinematic bifurcation of the MS. The observed 
morphology and kinematics suggest that three of these streams could be interpreted as a 3-way splitting 
of the main MS filament, while the fourth stream appears much younger and may have originated from 
the Magellanic Bridge. 

We find an extensive population of HI clouds at the tip of the MS. Two thirds of clouds have an angular 
size in the range 3.5'-10'. We interpret this as being due to thermal instability, which would affect a 
warm tail of gas trailing through the Galactic halo over a characteristic timescale of a few Myrs to a few 
hundred Myrs. We show that thermal fragments can survive in the hot halo for a long time, especially if 
surrounded by a < 10® K halo gas. If the observed clumpy structure is mainly due to thermal instability, 
then the tip of the MS is at a distance of ~ 70 kpc. A significant fraction of HI clouds at the tip of the 
MS show multi-phase velocity profiles, indicating the co-existence of cooler and warmer gas. 

Subject headings: ISM: clouds — ISM: structure — radio lines: ISM 



1. INTRODUCTION 

The Magellanic Stream (MS), a 10° wide tail of neutral 
hydrogen (HI) emanating from the Magellanic Clouds and 
trailing for almost 100° o n the s ky (Dec ^ -60 to +2Q°) 
(jWannier fc WrixonI [l972. : Mathewson et al.lll974l ). is the 
only clear example of a gaseous halo stream in the Milky 
Way's close proximity. While it is well accepted that the 
MS is the result of interactions between the Milky Way 
(MW) and the Magellanic Clouds, the relative importance 
of tidal stripping and various kinds of gasdynamical inter- 
actions is still very much under debate. Most recently 
several new attempts were made to model the observed 
HI column density and v elocity distribution as being due 
to purely tidal stripping (jConnors et alj 20061) . or gravita- 
tiona l -I- hydrodynamical interactions ( Mastropietro et all 
l2005f ). These models focus on reproducing general features 
in the Magellanic System, and gradients in HI column den- 
sity and velocity along the MS. To add more excitement 
to this topic, recent estimates of the proper motion of 
the Sm all Magellanic Cloud (SMC) bv iKallivavalil et al.l 
(|2006l ) and the m ost recent ca l culati ons of the Magellanic 
Cloud's orbits bv iBesla et al.l (|2007( ) suggest, contrary to 
all previous studies, that the Clouds are only on their 
first passage around the MW. The new orbits imply that 
neither tidal nor ram pressure stripping would have had 
enough time to produce the MS, calling for alternative for- 
mation mechanisms. 

The distance to the MS, especially to its tip, or the 
region the farthest away from the Magellanic Clouds, is 



another contentious question and varies greatly between 
models. Under the ram pressure hypothesis the tip has 
fallen the farthest toward the MW and is at a distance 
of on l y 25 kpc. Early tidal niodels ([Gardiner fc Noguchil 
119961 : lYoshizawa fc Noguchil Il999l ) place the tip at a 
distance of 60-70 k pc. The latest tidal simulations 
([Connors et al.ll2006f ) find an even more distant component 
exten ding from 170 to 2 00 kpc. The latest orbit calcula- 
tions ([Besla et al.l 1200*^ would also imply a significantly 
large distance to the MS, ^ 150 kpc, although it is not 
clear where exactly the MS is relative to the Clouds in 
this framework. 

For many years the MS was viewed as a complex of six 
discrete concentra tions (labeled as MS I to VI). New HI 
Parke s surveys by iPutman et al.l ([2003D and iBriins et al.l 
((2005'), with an angular resolution of 15.5', revealed a 
more complex nature of the MS gas, with a fascinating 
network of filaments and clumps. Two large spatial fila- 
ments were found to run in parallel over most of the MS 
length. Around Dec ~ 0° the dual filaments disperse 
into many small clumps and filaments culminating in a 
chaotic appearance at the tip. The only high resolution 
view of two selected regions at the tip of the MS was by 
Stanimirovic et al. (2002) who used the Arecibo telescope 
to image two small regions in MS V (Dec ~ 8°) and MS 
VI (Dec 12°). This work showed that the MS clumps 
have a complex morphology at 3.5' resolution, strongly 
suggestive of interactions between the MS and an external 
medium. While most previous studies thought that the 
MS dissipates at its tip (Dec ~ 0°), very sensitive West- 
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erbork observations by iBraun fc Thilkeil (|2004l ) suggested 
that the MS remarkably extends further to the north all 
the way to Dec 40°. 

Another interesting phenomenon brought to light by 
the Parkes surveys is the presence of numerous small HI 
clumps which sur round the main MS filaments in posi- 
tion and velocity ([Putman et alj [20031 ). While the origin 
of these clumps is still unclear, several possible mecha- 
nisms have been invoked: the dumpiness in the original 
gas drawn out of the Magellanic Clouds, the instabilities 
along the MS's edge, or dense condensations within an ex- 
tended mainly ionized MS component. One of the crucial 
issues about the origin and structure of the MS in general, 
is to what extent interactions with the MW halo determine 
or influence the MS gas. This problem becomes particu- 
larly important at the extreme northern end of the MS, 
because this portion of the MS is considered to be the old- 
est and has been immersed in the hot MW halo for a long 
time. 

In this paper we present results from the recent HI ob- 
servations of the tip of the MS obtained with the Arecibo 
telescope^ as part of the on-going survey by consortium 
for Galactic studies with the Arecibo L-band Feed Array 
(GALFA). In Section 2 we briefly outline our observing 
and data processing methods. Section 3 describes several 
new filaments discovered at the tip of the MS, while in 
Section 4 we present a statistical summary of properties 
of numerous HI clouds found in this region. In Section 5 
we discuss the origin of the MS and its clumpy structure 
based on our observations and results. 

2. OBSERVATIONS AND DATA PROCESSING 

The observations were conducted with the Arecibo tele- 
scope. GALFA HI survey consists of many individual 
projects. The data presented in this paper represent a 
combination of three GALFA projects: TOGS (or 'Turn- 
On GALFA Spectrometer' being undertaken in parallel 
with the ALFALFA extragalactic survey, PI: Putman), 
a2172 (Pis: Heiles & Peek, whose original target was a 
shell-like structure at Galactic velocities), and a2032 (PI: 
Stanimirovic, which observed a small region at the MS 
tip). 

As general GALFA o bserving and data reductio n strate- 
gies are summarized in [Stanimirovic "eTall (l200l . we em- 
phasize here only a few important points. All obser- 
vations were obtained with the dedicated spectrometer 
GALSPECT which has a fixed velocity resolution of 0.18 
km s~^. For a2032 and a2172, the telescope was driven 
in the basket-weave mode producing inter-woven zig-zag 
scans in the RA-Dec coordinate frame. For TOGS, nu- 
merous drift scans along right ascension were made, as 
this is t he preferred observing strategy for the ALFALFA 
survey (jGiovanelli et al.ll2005f ) . Before each observing scan 
about 5 minutes were spent running "The L east-Squares 
Frequency Switching" calibration procedure (|Heile£ll2007[ ) 
which allows us to derive a reference spectrum to be ap- 
plied to all spectra within the given scan. 

We have then combined scans from TOGS (the 'Fall' 
session) and a2172 to obtain the data cube that covers a 
region 37.5° x 24° in size. Various zig-zag and drift scans 

^ The Arecibo Observatory is part of the National Astronomy and 
agreement with the National Science Foundation. 



cross in many places. The crossing points were used to 
compare spectra and fine-tune our calibration. This pro- 
cedure removes small gain offsets between adjucent drift 
scans and produces a resultant image mostly free of strip- 
ping artifacts. After this fine-tuning, spectra from the two 
projects were gridded into a data cube. In the case of 
a2032 observations our strategy was different. After basic 
calibration, we have produced an a2032 data cube, and 
then added it linearly (in the image domain) to the com- 
bined TOGS-|-a2172 cube. Both combination processes 
worked very well. The images presented in this paper are 
the first demonstration of how successfully we can piece 
together data from different GALFA projects. 

The only problem encountered during the data reduc- 
tion is radio interference present is several TOGS scans 
at high negative LSR velocities. This interference affects 
all ALFA beams, varies across a given scan, and is cen- 
tered around 1422 MHz. It is caused by a still unidentified 
source. We have edited affected scans where possible, how- 
ever some artifacts are still noticeable in several images. 
Our future software developments will be better able to 
deal with this, and similar problems. 

As the integration time was different for the three ob- 
serving projects, our sensitivity varies across the image. 
ALFALFA and TOGS cover the sky twice with about 15 
sec of integration time per beam, while a2172 and a2032 
observations were obtained in a fast-scanning mode with 
an integration time of only about 5 sec per beam. Dec- 
hnation bands from 11° to 16° and 24° to 28° have the 
lowest rms noise level, 0.03 K per 1.4 km s^^ wide chan- 
nels. The rest of the image is noisier, with about 0.06-0.07 
K per 1.4 km s~^ wide channels. Our best 3-a sensitivity 
is 3.3 X 10^* cni"2 Q^g^. 20 km s"^ 

3. FOUR FILAMENTS AT THE TIP OF THE MS 

Figure [T] shows HI images from the final data cube at 
several selected velocities. In Figure [5] we show the first 
moment image, or the intensity-weighted velocity along 
the line-of-sight. The area covered in our observations 
is huge — about 870 square degrees (65°< I < 105°, 
-54°< b < -20°), focusing on the very tip of the MS 
where the two traditional MS filaments were found to be 
significantly spatially bifurcated (one is located at RA 23^ 
30"", and the other one at RA -23^ 00™). The only pre- 
vious high resolution observations of about 10 square de- 
grees within this region were obtained by Stanimirovic et 
al. (2002). 

As Figures [Hand [2] show, several coherent HI filaments 
extend all the way to Dec 25°. Braun & Thilker (2 001) 
reported an extension of the two MS filaments, in the 
form of diffuse gas, from Dec = 20° to 40°. We con- 
firm the continuous extension of the MS from Dec = 5° to 
~25°. Beyond this point our observations may not be sen- 
sitive enough to detect a very low-co lumn-density emis- 
sion found bv iBraun fc Thilkeil (|2004D . Our sensitivity of 
3 X 10^^ cm~^ (over 20 km s~^) is significantly worse than 
[Braun fc Thilkej (|2004[ )'s lowest contour (in their Figure 
5) shown at 3 X 10^^ cm~^. However, the picture of the two 
interwoven filaments changes dramatically in our observa- 
tions: instead of two traditional filaments we find several! 

Ionosphere Center, operated by Cornell University under a cooperative 
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Fig. 1. — Right ascension-declination images of HI emission at the tip of the MS for differentLSR velocities, given in the top left corner of 
each panel. The gray-scale intensity range is 0.005 to 0.15 K, with a linear transfer function. Several horizontal strips around Dec = 6°and 
13° , especially noticeable at lower negative velocities, are due to interference. The gap between Dec ~ 8° and 10° is due to lack of data in 
this region. The zig-zag pattern noticable in regions with a higher noise level is due to our scanning technique. 
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Fig. 2. — First moment image, or the velocity field, of the MS tip. Color represents the velocity centroids and brightness represents 
integrated intensity. The image has been smoothed to emphasize weak features. 





Fig. 3. — (left) An HI integrated intensity image with schematic representation of the four filaments, (right) Second moment or velocity 
dispersion. 
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While filaments often break into numerous branches and 
small clouds when viewed with Arecibo's resolution, there 
are at least four prominent and coherent large-scale struc- 
tures noticeable across the whole velocity range. To facili- 
tate further discussion we define and label filaments in the 
following way. Starting from the lowest negative velocities 
( 320 km s^i) in Figurc[Il 

f . 'Sf ' is centered at RA 23*^ 35™, this is a portion of 
one of the two traditional MS filaments; 

2. '82' is centered at RA 23'' fO" and can be traced 
through the whole velocity range; 

3. 'S3' is centered at RA 23'' 05"' and starts at a higher 
negative velocity of —360 km s^^; 

4. 'S4' starts at RA 23'' 00"' and extends to RA 22'' 
35"". 

Figure [3] (left) shows a schematic representation of the 
four streams. Both SI and S4 can be seen from the lowest 
negative velocities up to about —370 km s~^. S2, S3 and 
S4 seem to have spatially a common point of origin, and 
are most likely related to the second traditional MS fila- 
ment. The morphology of the four streams is significantly 
different: S2, S3 and S4 consist of and are surrounded 
by numerous small, compact clumps, having a beads-on- 
string appearance. The number of compact clumps peaks 
at velocities —350 to —370 km s~^. SI on the other hand 
has a more diffuse morphology and contains the smallest 
number of compact clumps. This may indicate that Sf 
has a different origin relative to the other three streams, 
or that it is younger and therefore did not have enough 
time to fragment. We get back to this issue in Section [51 

While generally clumpy morphology is prominent, the 
four defined streams are coherent structures, which can 
be seen most clearly in the first moment image. Figure [2] 
All four streams show clear velocity gradients, although 
not all gradients are similarly steep. In Figure [3] we show 
position-velocity images along the streams. SI and S3 
have a velocity gradient of about 5 km s~^deg~^ over a 
length of 15 degrees, in agreement with the overall veloc - 
ity gradient along most of the MS (' Putman et al.l [20031 ). 
S2 has the steepest gradient of almost 10 km s~^deg~^ 
along the first 10 degrees, and then it almost turns over 
in velocity roughly staying around —380 to —400 km s^^. 
S4 has almost a constant velocity of ^ —340 km s^^. All 
streams have a similar velocity dispersion, 1 5-20 km s~^, 
see Fi gure [3] (right). This is in agreement with lBriins et al.l 
(|2005[ ) who measured a dispersion of 15 km s~^ along most 
of the MS. 

4. COMPACT HI CLOUDS AT THE TIP OF THE MS 

A striking result from the above figures is the pres- 
ence of numerous small HI clouds. The smallest clouds 
have a size down to our resolution limit of 3.5' and are 
unresolved. Most of the clouds are compact and round, 
head-t ail morphology is su rprisingly uncommon. Previ- 
ously, iPutman et al. (|2003') noticed two types of clouds 
which surround the MS along its southern portion: head- 
tail clouds with tails pointing along the major axis of the 
MS (typically, the head has a column density of 10^^ cm~^ 
while the tail is 3-4 times fainter), and ubiquitous dense 



HI clouds which follow the main MS filaments in position 
and velocity, some again with a head-tail morphology. 

4.1. Cloud catalog 

We have compiled a catalog of 180 HI clouds. The po- 
sitions of all clouds are shown in Figure [5] Clouds were 
selected by eye from the spectral line data cube. The in- 
tegrated intensity image was used as a guide for cloud 
selection, local peaks in this image were considered as po- 
tential cloud candidates. We then used the full data cube 
to follow potential clouds spatially and in velocity, and 
measure cloud properties. Special care was taken to avoid 
breaking larger clouds into separate entities, however the 
catalog may still be slightly biased toward smaller clouds. 
We have made several checks to ensure cloud reliability 
and to exclude clouds potentially affected by scanning arti- 
facts or RFI. There were several cases (mainly at low Dec) 
where multiple clouds are present along the line-of-sight 
at clearly different velocities. These cases were treated as 
individual clouds. 

For each cloud we have measured the cloud size, the cen- 
tral velocity, and the peak HI column density. In addition, 
at the center of each cloud we have fitted the HI velocity 
profile with one or two Gaussian functions. Although the 
cloud catalog is not 100% complete, as it is hard to select 
clouds in the regions with higher noise, it certainly repre- 
sents well the cloud population whose physical properties 
we want to study. Basic cloud properties are presented 
and discussed in the following sub-sections. 

To ensure that we are dealing only with clouds of Mag- 
ellanic origin, our cloud catalog was cross-correlat e d wit h 
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the catalog of HVCs by iWakker fc van WoerdenI 
the catalog of HVCs and CHVCs bv lde Heii et all 
and the catalo g of recently discovered mini-HVCs by 
iHoffman et aP (I2004D . HVCs come with a wide range 
of sizes and properties. While some represent large com- 
plexes, like the MS, most are smaller with sizes of a few 
tens of degrees. On the other hand, CHVCs represent a 
class of HVCs which are compact (with sharp edges in the 
HI column density distribution) , apparently isolated from 
larger complexes, and with angular size less that 2°. One 
interpretation is that CHVCs are essentially the same as 
HVCs but located at much greater distances. Recently 
discovered mini-HVCs are even smaller and have diame- 
ters of 9' to 35'. These clouds are also more diffuse than 
typical HVCs, having the peak HI column density of only 
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Several of mini-HVCs were found to be 



superimposed on the northern extension (Dec > 0°) of the 
MS but at velocities distinct from those of the MS gas. 
One suggestion i s that these clouds r epresent outliers of 
the MS complex (iHoffman et al.ll2004D. 

There are 8 HVCs from Wakker fc van WoerdenI (|1991h 
in our field of view. Three HVCs are clearly de- 
tected and identified as clouds in our catalog, these are 
WvW413, WvW491 and WvW485. Additional three 
HVCs (WvW508, WvW492 and WvW472) have clouds 
nearby but at an angular separation larger than the beam 
of the Leiden-D wingeloo surv e y. Th ere are 19 HVCs and 
3 CHVCs from Ide Heii et all (|2003) in our field of view. 
All of these CHVCs were detected. 

As CHVCs were classified as isolated clouds it will be 
surprising to find a large population of CHVCs in our 
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Fig. 5. — A peak brightness temperature image of the whole region with cataloged clouds. Clouds without multi-phase structure are shown 
as crosses, while clouds with multi-phase structure are shown as circles. See Section 4.2. 
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Fig. 6. — Histograms of measured cloud properties: cloud angular size (FWHM in arcmin), the peak HI column density (in 10^® cm"'^), 
and the central LSR velocity. The dashed line in the left panel shows the angular resolution limit, while the dashed line in the middle panel 
shows the S-a sensitivity limit of the survey. 
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Fig. 7. — (left) Number of clouds as a function of Galactic latitude divided by the angular area observed at each latitude, (middle) Cloud 
size as a function of Galactic latitude. Dashed line represents a median value for each latitude bin. (right) Cloud central velocity as a function 
of Galactic latitude. In the case of clouds with multi-phase structure, the mean velocity is plotted. 
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cloud catalog. Only 3 CHVCs from Ide Heii et all (|2002[ ) 
are in our field of view. CHVC229 (from de Heij et al. 
2002) is extended and was found to be associated with sev- 
eral of our clouds at very low negative velocities (~ —330 
km s^^). We have excluded those HI clouds from our fur- 
ther analyses as at ~ —330 km s~^ these clouds are clearly 
spatially located far away from the MS emission. 

Three mini-HVCs are in our field of view: LMPU 
387-230, LMPU 387-275, and LMPU 387-297. LMPU 
387-275 and LMPU 387-297 are at velocities (-275 and 
—297 km s~^, respectively) where we detect several small 
clumps, the closest one being about 24' aw ay. We did not 
detect LMPU 387-230 however, although iHoffman et"ari 
(|2004 ) estimated its HI column density to be 4.8 x 10^^ 
cm~^, clearly above our detection limit. 

In summary, we do not find that our catalog is contam- 
inated by a significant population of CHVCs, it therefore 
contains primarily clouds associated with the MS. 

4.2. Cloud properties 

4.2.1. Cloud size, HI column density, and HI mass 

Figure [6] shows histograms of cloud basic properties: an- 
gular size, HI column density, and central (LSR) veloc- 
ity. The observed cloud angular size distribution strongly 
peaks around 10'. 90% of clouds have a size in the 
range 3' to 35'. The abundance of very small, compact- 
looking clouds is significant, more than half of clouds in 
the catalog have an angular size less than 10'. Obviously, 
these clouds would have been missed in previous large- 
scale su rveys such as the Le i den- Argentine-Bonn (LAB) 
survev (iKalberla et"an 12 005: Hartm ann fc BurtonI Il997t 
iBaiaia et al.l 20051 : Arnal et al..,2000,l . or even Parkes sur- 
veys of the MS. The low-size cut-off of the histogram is 
clearly affected by our resolution limit of 3.5', and higher 
resolution observations will probably find many smaller 
clouds. 

The cloud central HI column density distribution peaks 
around 1.3 x lO^^ cm'^, 90% range is 2.5 x 10^^ to 3.2 x lO^^ 
cm~^. This distribution function is affected by our sensi- 
tivity limit, deeper observations will probably find many 
clouds with lower HI column density. For comparison, 
the HI column density distribution of mini-HVCs peaks 
at 3 X lO^** cm-2. 

Figure[7](left) shows that the number of clouds decreases 
steeply along the MS, in the direction toward the tip. The 
number of clouds in each latitude bin has been divided by 
the observed angular area to account for different areal 
coverage across the image. Figure [7] (middle) shows a 
tendency for clouds with angular size larger than 20' to 
be found closer to the Magellanic Clouds (at b < —40°), 
while clouds with angular size less than 20' appear along 
the whole latitude range. This could be interpreted as ev- 
idence that part of the MS with b < —40° is at a larger 
distance, in agreement with tidal models. Please note that 
the gap around b = —25° on this plot is due to the lack of 
clouds, not due to the lack of data in this region. 

To provide an example of cloud linear size, HI mass, 
volume density, and thermal pressure, we use an angular 
diameter of 10' and a median peak HI column density of 
1.3 X 10"'^^ cm~^, for three representative distance values, 
20, 60, and 150 kpc, see Table [TJ 



4.2.2. Cloud central velocity 

Interestingly, the cloud central velocity, shown in Fig- 
ure [6] (right), shows a clear dichotomy. There are two 
peaks centered at —406 and —351 km s^^, with the ve- 
locity FWHM of 10 and 45 km s^^, respectively. As 
shown in Figure [7] (right) , the velocity dichotomy mani- 
fests in the way that there are two main groups of clouds, 
one around b = -47° to -37° (or Dec = 10°-15°) and 
V = —350 km s~^, and the other one around b = —35° to 
-27° (Dec = 18°-25°) and v = -400 km s-\ This figure 
also shows an almost continuous linear decrease in velocity 
(from -300 to -420 km s'^) with Galactic latitude. This 
is in agreement with the general velocity trend along the 
MS. 

The observed velocity dichotomy may indicate a kine- 
matic bifurcation of the MS tip. Such bif urcation would 
not be surprising. N-body simulations bv IConnors et al.l 
(|2006( ) found two kinematically separate components, one 
that follows closely the orbit of the LMC and has a lower 
negative velocity, and the other one with a more negative 
velocity. While separated kinematically, both components 
occupy the same (/, b) range. 

4.2.3. Clouds with multi-phase structure 

About 20% of the clouds in our catalog have veloc- 
ity profiles whose fitting requires two Gaussian functions. 
This is indicative of the existence of a multi-phase medium. 
However, some of these velocity profiles could be due to 
the existence of multiple clouds along the line of sight, 
whose velocity profiles partially overlap. To account for 
this effect, we accept as multi-phase clouds only those 
which have a subsonicly moving core and a warmer en- 
velope. This criterion translates into the requirement that 
the absolute difference between the centers of the two fit- 
ted Gaussian functions docs not exceed 0.7x the FWHM 
of the broad component. Our final sample of multi- phase 
clouds contains 21 entries. As an example, in Figure [S] 
we show velocity profiles of ten multi-phase (the first two 
rows), and five clouds with a single Gaussian profile (the 
last row). In most multi-phase cases, a narrow Gaussian 
component fits the line center, while a broader Gaussian 
function is needed to fit the line wings. 

Narrow line pro files were sca r cely fo und in the MS in 
previous studies. iBriins et al.l ()2005D found one cloud 
with a veloc ity FWHM of onl y 4 km s-\ at I = 93.6°, 
b = -53.7°. IKalberla fc HaudI (I2006D used the LAB sur- 
vey and decomposed all HI velocity profiles into Gaussian 
components. They found that the portion of the MS close 
to the Magellanic Clouds (and with positive velocities) 
has 27% of the sight lines with a narrow velocity compo- 
nent, while the portion further away from the Clouds (with 
negative velocities) has a significantly lower fraction, only 
about 10%. The lower number of narrow components was 
attributed to either a lower halo pressure and/or a larger 
distance of the northern MS extension. Our final frac- 
tion of multi-phase clouds (12%) is similar to Kalberla & 
Hand's estimate for the negative velocity portion of the 
MS. 

Figure [H] (left) shows the FWHM vs central velocity of 
all fitted Gaussian functions: crosses correspond to clouds 
with a single Gaussian component, while diamonds and 
circles represent the secondary (or narrow) and the pri- 
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Table 1 

Cloud radius, HI mass, HI volume density, and thermal pressure. 



Distance 


Radius 


HI mass 


HI volume density 


Pressure 


(kpc) 


(pc) 


(Mo) 


(cm~^) 


(K cm-3) 


20 


30 


2 X 10^ 


0.4 


40 


60 


90 


2 X 10^ 


0.1 


10 


150 


220 


1 X 10* 


0.05 


5 
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Fig. 8. — A few examples of HI clouds with and without multiple velocity components. The top two rows show clouds whose velocity 
profiles require two Gaussian functions, the bottom row shows a few clouds with just a single Gaussian velocity profile. For multiple velocity 
components, the broad velocity component is shown in green, the narrow component is shown in blue, and their sum is shown in red. 
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mary (or broad) velocity component for clouds with multi- 
phase medium. The FWHM of single-component clouds is 
predominantly in the range of 20-30 km s~^. For multi- 
phase clouds, the FWHM of the narrow component has a 
broad, almost continuous, range from 3 to 20 km s^^, with 
a median value of 13 km s~^. This would imply an upper 
limit on the kinetic temperature in the range 200-9000 K, 
mostly being in the t hermally unstable regime. For com- 
parison, iKalberla k. Hand ( 2006) find the FWHM of the 
narrow component in the range 5-28 km s^^. The FWHM 
of the broad component of multi-phase clouds is around 25 
km s""'^, in agreement with what is found for clouds with 
only a single Gaussian component. 

The right panel in Figure [9] shows the central velocities 
for both types of clouds. The velocity dichotomy is still 
noticeable with two peaks around —405 and —350 km s^^. 
This shows that the velocity bifurcation affects both clouds 
with and without multi-phase structure. We return to this 
point in Section r5. 2. 41 

5. DISCUSSION 

5.1. Origin of the MS filaments 

There are several important observational points pre- 
sented in the previous sections: 

• there are four coherent streams at the tip of the 
MS; 

• SI is diffuse, while S2, S3 and S4 are clumpy; 

• spatially, S2, S3 and S4 appear to originate from a 
similar location; 

• the velocity gradient along streams gradually de- 
creases from S2 to S4; 

• there is a lack of large clouds at low Galactic lati- 
tudes; 

• central velocities of clouds with and without multi- 
phase structure show velocity bifurcation at —350 
and —405 km s^^. 

Traditionally, it was as sumed that the MS originated 
purely from the SMC gas. iPutman et al.l (|2003t ) suggested 
that the two double-helix-like MS filaments could arise 
from gas stripped from the Magellanic Bridge and the 
SMC. By using th e LAB survey of Ga lactic HI (angular 
resolution of 0.6°). lNidever et al.l (|2007[ ) were able to trace 
one of the two filaments all the way to an over-density re- 
gion in the LMC. 

The only numeric al study that addressed spatial bifurca- 
tion of the MS is by[C onnors et all ()2006[ ). These authors 
showed that gross features of the MS and the leading arm, 
as well as the spatial bifurcation of the main MS filament, 
can be reproduced purely by tidal interactions between 
the MW and the Magellanic Clouds. They also suggested 
that the main MS filament is broken into two kinematic 
components. Several distinct features are found in their 
simulation. 

1. The main MS filament was created 1.5 Gyr ago in a 
major encounter between the SMC, the LMC, and 
the MW. Further tidal 'kicks' from encounters with 



the LMC 1.05 and 0.55 Gyrs ago resulted in its spa- 
tial, and then kinematic bifurcation. While one of 
the bifurcated components follows the orbit of the 
Magellanic Clouds and has a lower negative veloc- 
ity (~ —350 km s~^), the other possesses a higher 
negative velocity (—500 to —400 km s~^). 

2. In addition to the main stream, which extends in 
distance from 60 to ^ 100 kpc, there is the third 
stream right at the tip which spatially coincides 
with the main stream but has distinct kinematics 
and is significantly further away (distance of 170 to 
220 kpc). This component was created in an en- 
counter between the SMC and the LMC 2.2 Gyrs 
ago, and is the oldest and the most distant feature 
in the Magellanic System. 

3. Two tidal tails were drawn < 200 Myr ago from 
the Magellanic Bridge and follow the main filament 
along most of its length. 

By drawing a parallel with this simulation, streams S2, 

53 and S4 could represent a 3-way splitting of the main 
MS filament. As we show in the next section, this gas has 
had enough time to cool and fragment and is therefore 
very clumpy. As suggested by the simulation, most of the 
spatial 3-way splitting could have occurred in the more 
recent past, 1 Gyr to 500 Myr ago, during close encoun- 
ters between the MS and the LMC. The fact that S2 to 

54 show a decreasing steepness in their velocity gradients 
may be suggesting that these filaments were separated in 
subsequent encounters. SI on the other hand, may repre- 
sent one of the more recently formed tidal tails, drawn out 
of the Magellanic Bridge < 200 Myrs ago. The fact that 
it is significantly less clumpy than the other three streams 
suggests that SI may be significantly younger, and have 
not had enough time to cool and fragment. 

The lack of large HI clumps at the higher Galactic lat- 
itude suggests that the far tip of the MS may be more 
distant. However, the fact that we find clumps with multi- 
phase structure almost along the whole latitude range im- 
plies that the far tip is probably not too distant. As dis- 
cussed in the next section, it would be hard to justify the 
existence of multi-phase medium at a distance of ~ 200 
kpc. Another disagreement with the simulation concerns 
the observed LSR velocity at the far tip, we find ^ —420 
km s~^, while the simulation predicts velocities down to 
-500 km s~^ 

However, the observed large-scale filaments at the tip of 
the MS, as well as the evidence for kinematic bifurcation, 
give strong support to the tidal formation scenario for the 
MS. While the observations are only qualitatively in agree- 
ment with the simulations, and many details still have to 
be worked out, it is extremely encouraging to see that 
detailed spatial and velocity structure of the MS can be 
reproduced purely by tidal interactions. Better observed 
coverage of this region, as well as high resolution observa- 
tions of the southern MS extension, will certainly help to 
constrain this and other possibilities in the future. 

5.2. Origin of Compact clouds 

At certain velocities the MS has a very clumpy mor- 
phology. One of the obvious questions is: what physi- 
cal processes are responsible for such a clumpy structure? 
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Fig. 9. — (left) Velocity FWHM vs central velocity for all fitted Gaussian functions. Crosses correspond to clouds with only one Gaussian 
component. Diamonds and circles show narrow and broad velocity component, respectively, in the case of clouds with multi-phase structure, 
(right) Histogram of central velocities for clouds with a single Gaussian function (dashed line), and clouds with two Gaussian functions (solid 
line). Central velocities of both Gaussian components were used to derive the solid-line histogram. 



In trying to understand the cloud origin we look briefly 
into three important physical processes that would affect 
a tidal stream of warm gas moving through a hot ambient 
medium. All three processes have their characteristic time 
and spatial scales. 

5.2.1. Thermal instability 

Thermal instability arises due to gas cooling and results 
in spatial fragments of the WNM, whose size is defined 
by the cooling length Acooi- This process operates on the 
cooling timescale, defined by TcooI = fcTW/(Anw), with A 
being the cooling function, k being the Boltzmann con- 
stant, and Tw and riw being the temperature and volume 
density of the WNM gas. If we consider a WNM gas that 
originated in the SMC, then = 5 x 10~^ cm""^ (based 
on the size and total mass of the SMC, also measurements 
in (jStanimirovic et al.l[l999t )). then for = 8000 K we 
estimate A 

cool 

100 — 200 pc and the cooling timescale of 
Tcooi ^ 20 — 30 Myr. For a lower value of n^, which is not 
unreasonable for the gas drawn out of the SMC outskirts, 
both Acooi a-iid TcooI will increase. Considering that most 
models suggest the MS age of 1-1.5 Gyr, thermal instabil- 
ity, with a typical timescale of a few tens to a few hundreds 
of Myrs, most likely have had a significant effect on the 
evolution of the MS gas. 

Rec ent numerical simulations by lAudit fc Hennebelld 
(|2005[ ) provide an example of thermal condensation of 
warm gas. These simulations focus on a collision of in- 
coming turbulent WNM flows and thermal condensation 
of WNM into colder gas. While not exactly applicable in 
our case, it is still very instructive to notice the proper- 
ties of typical thermal fragments. A collision of incoming 
WNM streams creates a thermally unstable region, which 
partially condenses into cold clouds. The thermally un- 
stable warm gas (with temperature of 200 - 5000 K) has 
a fllamentary morphology and its fragmentation into cold 
clouds (with temperature < 200 K) is heavily controlled 
by turbulence. Cold clouds formed in simulations are ther- 
mally stable and long-lived, in the case of stronger turbu- 
lence they have round morphology, while a weaker turbu- 
lence is responsible for more elongated morphology of cold 



clouds. The simulations show that the fraction of cold 
gas (with temperature < 200 K) produced in this process 
ranges from 10% in a strong turbulent case, to about 30% 
in a weak turbulent case. It is interesting to note that the 
lower bound on this CNM-to-WNM fraction agrees with 
the fraction of HI clouds with cool cores we flnd in our 
observations (12%). 

5.2.2. Kelvin- Helmholtz (KH) instability 

Kelvin-Helmholtz (KH) instability occurs at the inter- 
face of the moving warm stream of gas and the ambient 
medium and also results in the formation of small frag- 
ments with a power-law size distribution. The timescale 
on which this instability develops, tkh (given in Gyr), is 
given by: 



(jMurrav et al.lflQQ^ . where Rd would be the initial size 
(radius) of the MS gas pulled out of the Magellanic 
Clouds, Da = n^/rih corresponds to the density con- 
trast of the warm KH fragments and the ambient hot 
medium, and v is the s tream velocity. We use v = 220 
km s~^ ([Stanimirovic et al. 2002). assuming that the ve- 
locity of the original MS gas is equal to its present day 
velocity. We also assume that Rc\ ~ 10°, based on the size 
of the MS VI portion of the MS, and that this material was 
initially at a distance of 60 kpc. 

If rih = 10""* cm~'^, and riw = 5 x 10^^ cm~^, then 
Tkh ~ 600 Myr. As this timescale is almost comparable 
to the age of the MS (in the framework of tidal and ram- 
pressure models), we can conclude that the KH instability 
has had a negligible effect on the evolution of the MS gas. 

5.2.3. Cloud evaporation 

Once the warm stream has fragmented due to ther- 
mal/KH instability into smaller clouds, the third process 
becomes important — this is conductive heat transfer be- 
tween th e warm stream fragm ents and the hot ambient 
medium. iLin fc Murray! (|2000f ) show that the size of the 
smallest fragments is set by the cloud ability to radiate 
away the energy flux due to conduction. 
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The warm clouds reside within warmer phases, hot halo 
gas with Th ~ 10*' K and warm ionized gas with Th ^ 10^ 
K. Detections o f OVI a bsorption from gas associated with 
the MS .Sembach et al.1 (^003) give strong support for the 
existence of an ionized component around the MS with 
Th < 10^ K. It is expected that different phases are sep- 
arated by an interface region through which heat flows 
between the two phases. The interface is a surface of 
transition where, for small clouds, the WNM evaporates 
into the warmer medium. The classical evaporation the- 
ory (|McKee fc Cowielll977[ ) predicts a critical radius for 
clouds at which radiative losses balance the conductive 
heat input, and at which the cloud neither evaporates nor 
condenses. Clouds smaller than the critical radius evap- 
orate, while clouds larger than the critical radius accrete 
by condensation of the surrounding gas. 

In the case of an interface between the WNM and the 
HIM, the cloud critical radius (i?rad) times the density of 
the surrounding medium, rth, is 5 x 10"'^^ cm~^ for Th — 10^ 
K, or 5 X 10^^ cm-2 for Th = 10^ K (from Fig. 2 in McKee 
& Cowie 1977). Generally, the expected critical radius de- 
creases with increasing rih. However, for nh ~ 10^* cm^'^, 
the critical radius is large, i?rad > 200 pc. Clouds smaller 
than ~ 200 pc would then evaporate, while WNM clouds 
larger than ^ 200 pc would accrete the HIM by condensa- 
tion. Thermal fragments we considered in Section 5.3 with 
a size Acooi ~ 200 pc will be predominantly smaller than 
the critical radius and therefore undergoing evaporation. 

These small (~ 200 pc), evaporating clouds are sur- 
prisingly long-lived though. For mainly neutral clouds 
surrounded by the hot gas, the classical mass-loss rate 
is 2 X 10^^ gr yr~^ for Th = 1 0*^ K, or ~ 10^^ gr yr'^ 
for Th = 10^ K (|Slavinl l2007h . In the case of satu- 
rated evaporation the mass-loss rates are even lower. For 
riw = 5x10"^ cm~'^ the evaporation timescale is r^vp = 0.4 
Gyr for Th = 10*^ K, or Tevp = 80 Gyr for Th = 10^ K. But 
the timescales are in general long, approaching or being 
longer than 1 Gyr. We hence conclude that small frag- 
ments would be able to survive for a long time, at least a 
Gyr, most likely throughout the whole lifetime of the MS. 

5.2.4. Constraining cloud distances 

As thermal instability must have played a very impor- 
tant role in the evolution of the oldest MS gas, we suggest 
that this could be the dominant process responsible for 
the clumpy structure found in our observations. The im- 
portance of thermal instability is further supported by the 
large abundance of small HI clumps (see Section 4.2.1). 
The KH instability is important as well, however it has 
a significantly longer timescale. As we have shown, it is 
reasonable to expect that thermal and/or KH fragments 
can survive for a long time, especially if surrounded by the 
ambient gas with Th < 10^ K. 

It is interesting that cloud morphology in our observa- 
tions is dominated by round shapes. Elongated or head- 
tail clouds are surprisingly absent in our HI images. Head- 
tail morphology is commonly observed for CHVCs and 
is often seen as a signature of ram-pressure stripping, or 
other types of interactions between a cloud and the ambi- 
ent halo medium. We therefore conclude that ram pressure 
effects probably have a secondary role, with gravitational 
forces being the dominant gas removal process. 



If we assume that thermal fragmentation is the domi- 
nant process responsible for the clumpy structure, we can 
estimate a distance to the MS. If we assume that the me- 
dian value of the cloud angular size distribution (10') cor- 
responds to the linear size of 'typical' thermal fragments, 
Acooi ~ 200 pc (derived assuming riw = 5 x 10~^ cm~^), 
we estimate a distance of ^ 70 kpc. Please note that these 
numbers are only representative. The cloud angular size 
varies and is also affected by our resolution limit, the vol- 
ume density of the warm MS gas is also not a constant 
number. However, this exercise is very instructive and 
demonstrates that thermal fragments can reasonably well 
explain at least some small-scale structure in the MS. The 
MS distance of ~ 70 kpc is in agreement with the pre- 
dictions from tidal models. Even more impressively, this 
distance agrees well with the recent estimate of 75 kpc by 
Jin & Lyndcn-Bcll (|2007f ). who used a geometrodynami- 
cal model and assumed that the MS follows a planar orbit 
around the Galactic center. 

Also, the existence of multi-phase HI clouds can be used 
to place an upper limit on the MS distance. It is im- 
portant to stress that the existence of multi-phase clouds 
along the tip of the MS, and especially at a distance of 
~ 80 kpc, is extremely surprising from a theoretical point 
of view. Based on the consideration of cooli ng and heat- 
ing processes in an isothermal 10^ K halo, Wol fire et al.l 
(1995) predicted that multi-phase clouds, pressure con- 
fined by the hot halo, should exist in the MS only at dis- 
tances less than 20 kpc fro m the Galactic p lane. More re- 
cently, considering CHVCs [Sternberg et al.l (j2002) showed 
that dark matter can provide an additional confinement 
mechanism, allowing the existence of multi-phase struc- 
ture at distances <j 150 kpc. The multi-phase clouds at 
the tip of the MS therefore suggest that the tip is likely 
to be at a distance < 150 kpc. This argues against the 
existence of th e third, very d i stant stream found in the 
simulations by IConnors et al.l (|2006t ). at least within the 
area covered by our obs ervations. Also, as pointed out by 
iKalberla fc: HaudI ()2006D . the multi-phase structure at the 
MS distance suggests a lower confining halo gas pressure 
than what is expected in the theoretical framework. 

6. SUMMARY AND CONCLUSIONS 

We have obtained HI observations of a region, ~ 900 
square degrees large, at the tip of the MS. The major re- 
sults are as follows. 

1. We find four coherent HI streams at the tip of the 
MS, extending continuously from Dec 5° to 25°. 
The streams have different morphology and veloc- 
ity gradients. Three streams are very clumpy, while 
one stream is mainly made up of large diffuse struc- 
tures. 

2. The compa rison of observations with the tidal sim- 
ulations W lConnors et al.l (|2006[ ) suggests that the 
three clumpy streams could be a result of a 3- 
way tidal splitting of the main MS filament, in- 
duced by encounters between the LMC and the MS. 
The fourth stream is probably much younger, and 
could have been tidally drawn from the Magellanic 
Bridge in the more recent past, < 200 — 500 Myr. 
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This stream is more diffuse because it has not had 
enough time to cool and fragment. 

3. We find an extensive population of HI clouds at 

the tip of the MS. The smallest clouds have a size 
reaching down to our resolution limit of 3.5'. 70% 
of clouds have an angular size in the range 3.5'-10'. 
About 12% of the MS clouds show multi-phase ve- 
locity profiles, indicating the co-existence of cooler 
and warmer gas. Cloud central velocities show two 
peaks at —350 and —400 km 

4. We show that for warm gas trailing through the 
Galactic halo for the past Gyr or so, the effects of 
KH instability are relatively unimportant. How- 
ever, thermal instability, with a typical timescale of 
a few Myrs to a few hundred Myrs, is significantly 
important and will result in the fragmentation of 
gas and the formation of smaller clumps. The ob- 



served dumpiness at the tip of the MS could be 
fully or partially due to this mechanism. 

5. We show that thermal/KH fragments can survive in 
the hot halo for a long time (> 1 Gyr), especially if 
surrounded by a < 10^ K halo gas. 

6. If the observed angular size of HI clumps is due to 
thermal instability, then the tip of the MS is at a 
distance of ~ 70 kpc. 
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